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ABSTRACT 

We investigate the effects of the observed UV background radiation on galaxy forma- 
tion. Photoionization by UV radiation decreases the cooling rate of the gas in halos, so 
that objects with only large density contrasts can self-shield against the background radi- 
ation, thereby allowing the shielded neutral cores to cool and form stars. In the context 
of the CDM model, we use the criterion that self-shielding is essential for star formation 
to calculate the mass function of galaxies based on both Press-Schechter and the peaks 
formalisms. The ionizing UV radiation results in the inhibition of galaxy formation — 
we show the decrease in the number density of the galaxies quantitatively. We also find 
that the merging in general is made inefficient by the UV photons through photoioniza- 
tion of the gas in the bigger system into which small objects are incorporated. The latter 
means that, in a merging dominated region, where the number density at the low-mass 
end (Mb ;$ 10 10 Mq) is usually expected to decrease with time, the trend is reversed (the 
number of low-mass galaxies is more at lower redshifts z) due to the decreasing UV flux 
with time after z ~ 2. We further discuss the implication of our results for the number 
counts of galaxies and possible evolution of the luminosity function of the galaxies. 

Subject headings: formation, evolution - galaxies ; diffuse radiation - cosmology 



1. INTRODUCTION 

Formation of galaxies depends sensitively on the dissipative nature of the baryonic 
gas. While the halos are formed out of the underlying distribution of the dark matter, if 
there is any, the fate of the resulting galaxy depends much on the behavior of the gas in 



the halo. Cooling of the gas has long been known to be a decisive factor in determining 
the characteristic size and mass of the galaxies (Rees & Ostriker 1977; Silk 1977). Any 
deviation in the history of the universe that alters the state of the gas in halos, therefore, 
must result in some change in the property and distribution of galaxies. 

One motivation to study such effects is the fact that standard model of hierarchical 
structure formation predicts more low mass galaxies than observed. One then requires 
a negative feedback mechanism to halt or inhibit the formation of small galaxies (with 
luminosities less than the characteristic value ~ 10 10 h~ 2 L@), so that most of the mass 
of the universe is not locked away in such galaxies. Two mechanisms have been put forward 
to achieve this: heating by supernovae from the first bout of star formation and heating 
by the UV background radiation. Dekel and Silk (1986) (and recently elaborated upon 
by Lacey & Silk 1991) have considered the first mechanism in detail. They found that 
galaxies with velocity dispersion less than ~ 100 km s _1 are vulnerable to substantial 
mass loss due to supernovae driven winds. With this limiting value of velocity dispersion, 
and assuming that dwarf galaxies form out of la density fluctuations, they derived several 
scaling relations and compared them with observations. 

The possibility that UV photons can inhibit the collapse of small galaxies has been 
recognized by various authors (e.g. Dekel & Rees 1987). However, attempts have been 
made only recently to measure, and calculate from theoretical arguments, the intensity of 
the UV background radiation. While the precise measurements are still eagerly awaited, 
it seems almost certain that a background radiation with intensity (at the Lyman limit) 
J ~ 10 -22 — 10 -21 erg cm -2 s _1 sr _1 Hz -1 permeated the universe at redshifts 5 £ z £ 2 
(see, for example, Bechtold et al. 1987; Espey 1993). The UV photons are thought to have 
been produced by quasars and, possibly, young galaxies (Miralda-Escude & Ostriker 1990; 
Madau 1992). 

Since gas with T ^ 10 5 ' 5 K cools mainly by line cooling of H and He + atoms, the 
cooling efficiency is diminished in the presence of the ionizing photons. Efstathiou (1992) 
calculated the cooling function and, after comparing the cooling time with the Hubble 
time, argued that galaxies with velocity dispersion ^ 50 km s _1 will not be able to cool 
and collapse. However, it is not clear as to how much the UV background affects the 
mass function of the galaxies , and the evolution of the mass function with time (since 
the background radiation changes with time). In the hierarchical model of structure for- 
mation, density perturbations collapse to form halos at different redshifts depending on 
the amplitude of fluctuation. Halos with very large amplitude, i.e. large density contrast, 
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will be able to self-shield themselves against the UV background and can collapse even if 
the velocity dispersion is small. Therefore, without detail calculations, one cannot answer 
the question - to what extent does the UV background inhibit galaxy formation, say, at a 
certain mass range, if at all? 

We believe that such questions are very important, since suppression of galaxy forma- 
tion via photoionization has been invoked in various contexts. To explain the paradoxical 
observations of galaxy counts at faint magnitudes (e.g. Cowie, Songaila, & Hu 1991), 
Babul and Rees (1992) have argued that small mass objects of mass ~ 10 9 M collapse at 
z ~ 3 but do not form stars until z ~ 1 due to the inhibiting UV background radiation. 
After the starburst, galaxies with dispersion velocity less than the Dekel and Silk limit of 
v ~ 100 km s _1 are destroyed so that not many of the initial set of galaxies are observed 
today. The effect of the UV background radiation on the abundance of small mass cfalax- 
ies, therefore, are important to understand in order to make any model for the faint blue 
galaxies. 

In this paper, we discuss the implications of the UV background in more detail. Our 
principal goal is to make the estimates and guesses of the previous workers for the effect 
of UV radiation on galaxy formation more quantitative. We take advantage of the recent 
measurements of the UV background to calculate the exact changes, some of them hitherto 
unforeseen, caused by the UV photons in the mass function of galaxies. For simplicity, 
and because the power spectrum in cold dark matter (CDM) model has been calculated in 
more detail than in any other model, we will use the CDM spectrum, with Q — 1. For the 
Hubble constant, we use /iso = 1, where H — 50/iso km s _1 Mpc -1 . The structure of the 
paper is as follows: in §2 we calculate the cooling function, and discuss self-shielding of 
clouds against UV photons. In the density-temperature plane, we plot the loci of constant 
amplitude density perturbations and determine which objects can form a neutral core 
and thus cool. In other words, we determine the threshold u t f l (M) for the amplitude of 
fluctuation as a function of the UV background intensity. In §3 we use the Press-Schechter 
and the peak formalism of Bardeen et al. (1986) to calculate the mass function of galaxies 
with and without the UV background. We discuss the implications of our results in §4. 

2. UV BACKGROUND RADIATION AND EVOLUTION OF GAS CLOUDS 

2.1 UV background radiation 
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It has been evident since long that the intergalactic gas is highly ionized at high 
redshift (Gunn & Peterson 1965). Webb et al. (1992) have recently claimed to have 
measured an optical depth of r ~ 0.04 at z — 4. Jenkins & Ostriker (1991) put an upper 
limit of the optical depth due to neutral hydrogen in the intergalactic medium of r < 0.14 
at z — 4.2. While it is possible that the gas could be ionized due to processes other 
than photoionization, observations of Lyman-a clouds independently suggest the existence 
of a UV background radiation from the 'proximity effect'. The accurate intensity of the 
ionizing radiation is yet to be confirmed: the reported flux varies from (in the units of 
J_2i, where J = J_ 2 i x 10 -21 erg cm" 2 s" 1 sr" 1 Hz" 1 ) ~ 1 (1.8 < z < 3.5) (Bajtlik et 
al. 1988) to ~ 0.1 {z ~ 2) (Dobrzycki & Bechtold 1991; Espey 1993). On the theoretical 
front, Miralda-Escude and Ostriker (1990) suggested that quasars and young galaxies could 
give rise to a flux of J-21 ~ 1. Madau (1992) calculated the contribution of the quasars 
and concluded that quasars solely could account for the ionizing flux if J-21 ~ 0.1 at 
z ~ 2. His calculations suggest that the intensity is almost a constant for z £ 2 and 
decreases rapidly at lower redshifts. For the present day value of J-21, we will use the 
observations of the sharp edge of HI in spiral discs by Maloney (1992, 1993), who inferred 
that J_2i (2 = 0) = 0.033 x 10 ±0 7 . 

Below, we will assume a spectral index of the UV radiation a — 1 ( i. e., J u oc (^) -1 )- 
Since Black (1981) has calculated the atomic processes for a — 1 (which we will use to 
construct the cooling function), and since, the expected quasar spectrum also has a similar 
index (a ~ 0.5), our assumption makes the calculation simple. 

2.2 Cooling rate 

The cooling rate of a plasma at T ^ 10 5 ' 5 K with primordial abundance (X=0.76, 
Y=0.24) depends mostly on the fractional ionization of H and He + . In the absence of 
any ionizing flux, the cooling is dominated by line emissions from H and He + atoms. The 
cooling is therefore highly diminished when UV photons exist (which ionize these atoms) 
and the cooling is then dominated by recombination. The cooling rates in the presence of 
ionizing radiation with intensity J-21 — 1, 10 were calculated by Efstathiou (1992) using 
the processes tabulated by Black (1981). We have also calculated the cooling rate using 
the same processes. We find a slight discrepancy, by a factor < 1.3, with that shown 
in Efstathiou's (1992) paper and realize that such discrepancies abound in the literature 
(compare, e.g., the cooling rate in Efstathiou 1992 for no ionizing flux with that in Fall & 
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Rees 1985, both calculated for a primordial gas). We should, therefore, bear in mind the 
uncertainties when making any argument based on the cooling rate. 

The heating rate is also determined by using the processes in Black (1981). The 
equilibrium temperature of the plasma in the presence of ionizing flux is calculated by 
balancing the heating and the cooling rate. 



Efstathiou (1992) plotted the limiting curves where cooling time, t coo i, equals the 
Hubble time, tjj, at various redshifts to show that photoionization can be important as 
an inhibiting mechanism. We have instead decided to plot curves comparing cooling time 
and the collapse time which are independent of redshifts (compared to various plots for 
different redshifts in Fig. 2 of Efstathiou 1992). 

Fig. 1 is a plot in the particle density (n) and temperature (T) plane. Cooling time, 
t coo i is defined as (Efstathiou 1992) 



where -L(T, n) (erg cm 3 s _1 ) is the cooling function, fi is the mean molecular weight (= 0.59 
for a totally ionized gas) and — /^(l — Y)n is the number density of hydrogen atoms. 
The time scale of collapse of an object, tff is given by (Dekel & Silk 1986), 



where we have assumed that the gas mass comprises only 10% of the total mass. 

The t coo i — tff curves for J_2i = 0,0.1,1 are shown in Fig.l (when comparing with 
Fig. 5 of Dekel and Silk (1986), note that their n is actually n# in our definition). 

We have plotted the curves for constant amplitude of the density fluctuations in CDM 
model. The curve marked 'la' denotes the perturbations for which SM/M — 1 at a 
comoving radius of 16 Mpc. We have used the fact that, if Sq is defined as the linear 
density fluctuation extrapolated to the present time, then Sq — vA(M), where A(M) is the 
variance (an analytical approximation is given in White & Frenk 1991). Spherical top-hat 
model gives, 



2.3 Cooling time vs. collapse time 




3 1 kT 
2/i 2 (l -Y) 2 nL(T, ra)' 



(1) 



t ff = (6nGp)- 1/2 = 8.9 x 10 14 n 



- 1/2 s, 



(2) 




(3a) 



(36) 
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where m p is the proton mass and po is the present mean density of universe with h$o — 1 
and Q — 1. These relations allow us to draw the constant amplitude curves in the n — T 
diagram. One can scale the curves for other values of the biasing parameter, b (we have 
assumed b — 1 for Fig.l), using Sq oc vjb. 

The corresponding curve (separated by an arrow) denotes the gas clouds after a con- 
traction by a factor of _F -1 / 3 , or equivalently the increase of density by a factor of F, where 
F is the ratio of gas mass to the total mass (on the assumption of uniform density). The 
gas densities on this curve are comparable to the halo densities so that star formation can 
ensue (White & Rees 1978; Dekel & Silk 1986). We shall next require that the gas clouds, 
after a collapse to this stage, be self-shielded against UV photon to ensue cooling and star 
formation. 

The density for the above curve is calculated for a uniformly distributed gas in CDM 
halos, with the fraction of gas mass being 0.1. The short-dashed diagonal lines denote 
objects of a given mass (again with a gas mass fraction of 0.1). 

The equilibrium temperature T eq is then plotted as a function of the density for 
J_2i = 0.1, 1. It is easily seen that for the objects that can cool faster than their collapse 
time, T eq ranges between 10 4 - 10 5 K. 

2.4 Self- shielding 

Large objects can shield themselves against the UV background radiation. For pertur- 
bations with large density, the optical depth to the UV photons can reach unity at a certain 
radius and the gas inside may be shielded against any UV photons and can, therefore, cool 
fast. There exists a critical radius for which the rate of photoionization (oc r 2 , the surface 
area) balances the rate of recombination (oc r 3 , the volume). The critical radius, r c , is 
given by the balancing equation (assuming a spherical and uniform gas cloud), 

47T 

—r 3 c n e n p a{T) « (4vrr 2 ) 1.5 x 10 5 J_ 21 , (4) 

where cx.iT) is the recombination rate for hydrogen atoms, and Tig , Tip clXC the electron and 
proton densities. We use T eq in the above equation, assuming that objects with larger 
virial temperatures will cool down (t coo i is much shorter than the Hubble time for objects 
inside the region t coo i <tff) and those with lower virial temperatures will be heated up 
to T eq quickly. 

The sphere with the critical radius, r c , is essentially an 'inverted Stromgren sphere'. 
Objects bigger than this will have a core where no UV photons can reach, because the 
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surrounding mantle of ionized gas balances the photoionization with recombination. (The 
r c calculated here refers only to hydrogen; however, there would be a thin shell where both 
H and He + will be ionized, surrounding a shell where H is ionized but He + will be neutral. 
We neglect such details, since recombination due to H is dominant.) 

The virial temperature of the object with the critical radius is easily calculated (total 
mass being 10 times the gas mass), when they are assumed to be self-gravitating. We 
plot the curves for these critical objects in Fig.l for J-21 = 0.1, 1 (the dotted-and- dashed 
lines) and will refer to them as the 'self-shielding curves'. 

One can calculate the core mass for objects lying above the self-shielding curves, by 
calculating the thickness of the shell, Ar, where recombination and photoionization are 
balanced. One can easily show that 

^R-WR- W™* 1 *'-" )] 1 '', (5) 
L V n e tip a(l ) /J 

where R is the radius of the gas cloud, and n e , n p , a are expressed in cgs units. Fig. 2 
compares the core mass (after a dissipative contraction by a factor i* 1-1 / 3 ) with the total 
mass for different values of J-21 and v. The points where the core masses plunge to zero 
are essentially the points of intersection of constant a curves after a contraction with the 
self- shielding curves in Fig.l. 

Note that in deriving eq.(5) we have only considered photoionization. Collisional 
ionization would be important for clouds with T eq > 10 5 K and would lower the self- 
shielding curves at the high temperature region of the n — T diagram. We will neglect this 
effect in the following and only note here that it would make our results more conservative. 

2.5 Implications of the n — T diagram 

The cooling diagram Fig.l shows the features that have been discussed qualitatively by 
previous works, viz., the suppression of the dwarf galaxies. If the dwarf galaxies originate 
from typical la fluctuations, as is commonly believed (Dekel & Silk 1986; Babul & Rees 
1992), then their formation will be inhibited by the UV photons if J-21 ~ 10. Figure 1 
shows that for J-21 ~ 1, objects with v ~ 0.5 will not be able to self-shield and cool, 
after collapsing by a factor of i* 1-1 / 3 . With the decrease of the UV flux after z ~ 2, the 
self-shielding curves will be lowered in Fig.l, and objects that were earlier above the curves 
will now be able to form neutral cores and harbor star formation. 

^From Figs. 1 and 2, it is quite obvious that for a given amplitude of the density 
fluctuation, v, objects with larger masses (and therefore with smaller densities) are more 
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vulnerable to the UV background radiation. These objects remain ionized and the heating 
due to the UV photons keep the temperature at £ T eq and cannot cool. In other words, 
one can say that for a given mass M, the value of v must be bigger than some lower limit, 
Vth{M). This threshold value of v is shown in Fig. 4 as a function of the total mass and 
J-21 (see the next section for details). Armed with the knowledge of v t h(M) we can now 
calculate the effects of the UV background on the mass function of galaxies in the next 
section. 

3. GALAXY MASS FUNCTION IN THE UV BACKGROUND 

In the presence of photoionizing radiation, the number of galaxies which can cool 
before a given epoch is reduced. Since the effect depends on mass, as we have seen in the 
previous section, the shape of mass function of galaxies will inevitably be changed by the 
UV background. Blanchard et aZ.(1992) have also discussed, using the different model, the 
effects of gas heating on the galaxy mass function. In this section, we describe the basic 
formalism to derive the mass function and discuss the results. We will not touch upon the 
luminosity function of galaxies, which can be deduced only by introducing models of star 
formation (e.g. Lacey & Silk 1991; White & Frenk 1991, Lacey et al. 1993). Since the 
exact theory of star formation still eludes us, we will only derive the mass function and 
attempt some model-independent conclusions for the luminosity function. 

3.1 Basic equations 

We calculate the mass function using Press-Schechter formalism (Press & Schechter 
1974, hereafter PS) and the peaks formalism (Bardeen et al. 1986). In the PS formalism, 
mass function of galaxies at a fixed time is determined by earlier history of hierarchical 
mergings — i.e. the number of galaxies is reduced by their mergers into larger collapsed 
systems. In the peaks theory, on the contrary, the peaks which will eventually be within 
the larger peaks are counted as independent objects — i.e. no merging is assumed. 

PS formalism gives an analytical description for the number of galaxies at a given 
epoch which have not been incorporated into the collapsed larger system. It applies to the 
evolution of dissipationless dark halos, which are subject to hierarchical merging as seen 
in cosmological N-body simulations (e.g. Frenk et al. 1988). It is however not obvious as 
to whether baryonic cores, once formed, merge together as rapidly as the dark halos do. 
The small cross sections for core-core interaction would reduce the rate of merging, so that 
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one-to-one correspondence between dark halos and baryonic cores may not be established 
(see, e.g. Carlberg & Couchman 1989; Carlberg 1993). This is the reason why we also use 
the peaks formalism in which no merging is assumed to estimate the mass function. The 
actual number of galaxies at a given mass probably lies between the values predicted by 
these two formalisms. 

In CDM model, it is assumed that galaxies arise from the linear density fluctuation 
5 = Sp/ p with Gaussian statistics, where p is the mean density of Universe. The idea is 
that if the mean linear overdensity extrapolated to z — is larger than some critical 6 C , 
then the (spherical) perturbation will have eventually collapsed and virialized. We adopt 
S c — 1.69 as given by the spherical "top-hat" collapse model. 

In the standard PS formalism, the number of objects per comoving volume per unit 
mass is usually written as, 



where po is the present mean density of Universe; v is the amplitude of density perturbation 
for the onset of collapse in the units of rms density fluctuation, and is defined as 



where z c is the epoch of collapse. A(M) is the rms linear overdensity for the mass scale 
M at the present time. Equation (6) applies to galaxies which have not been incorporated 
into larger collapsed system at a given z c . Here we note that one must be careful about this 
standard form: eq.(6) is derived on the assumption of constant critical 6 C . However the 
adoption of the cooling and self-shielding criterions as we will consider below introduces 
the mass (or equivalently, the filter radius)-dependent critical overdensity 6 C (M), and the 
simple incorporation of this effect into eq.(6) is beyond its range of validity. This invalidity 
of the standard PS formula in general cases of 6 C (M) is related to the factor of 2 fudge to 
derive eq.(6): half the mass which is uncounted in PS, i.e. underdense regions with 5 < 6 C , 
can have 5 > S c for some larger filter in general. 

The explicit consideration of this effect in the PS formalism was given by Peacock & 
Heavens (1990) and Bond et al. (1991), and we will adopt the former recipe, which can 
apply to general cases of 6 C (M). Denoting by pg(£> > <*>c) the usual Gaussian probability 
that a point has 5 > S c for a given filter radius Rf, Peacock & Heavens showed that the 
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probability p(Rf) that the point is ultimately incorporated into a larger Rf system is given 

CIS j 

p(Rf)=PG+(l-PG)(l-Ps), (8) 

where p s (Rf) is the survival probability that the point is always below S c for larger Rf, 
given in eq.(ll) and (14) by Peacock & Heavens (1990). The differentiation of this by mass 
defines the mass function per unit volume 

Mn(M,z)dM = Po \^\. (9) 

On the other hand, in the peaks formalism, the differential number density of peaks 
for a given Rf is given as 

dripk _ 1 



dv (2vr) 2 i?3 



exp(-i/72)G( 7 , 7 i/). (10) 



where R* and 7 are defined in terms of a set of spectral moments a t , i.e. R* — \f^>o\jo2, 
and 7 = o\jo^02i respectively, and G is defined in Bardeen et al. (1986). For the differ- 
ential number density in the peak size Rf, we adopt the approximate approach given by 
Lacey & Silk (1991): 

d 2 n pk 1 3 / Rf\ 3 dlnR* 2 , , , , 

exp(-i/ /2)G(7,7i/), (11) 



dvdRj {2n) 2 R f \R*J dlnR f 

and integrate this over the peak height v to obtain the mass function. 

We use the power spectrum derived by Davis et al. (1985). It is important to notice 
that if we adopt the top-hat filter in these two formulations, which has a sharp edge in 
r-space but is oscillatory in &-space, the integrals for spectral moments do not necessarily 
converge for all possible slopes of the fluctuation spectrum. The Gaussian filter avoids 
this undesirable feature, and thus we adopt it, using the same power spectrum. In order 
to relate the ingredients of top-hat collapse model, such as S c — 1.69 for collapse, to the 
Gaussian filtering, we have normalized the Gaussian-smoothed A(M) to be unity at the 
corresponding top-hat lengthscale R of the usually adopted value 16/i^" 1 Mpc. The mass of 
a Gaussian-filtered peak is estimated by M — (27r) 3 / 2 R^po, so that the relation between 
R and R f is given by R = (3/47r) 1 / 3 (27r) 1 / 2 i? / . In Fig. 3, we have plotted the two A's, 
as functions of the top-hat lengthscale. The difference between the two cases is small. 
Of course, assigning a proper volume to the Gaussian filter function in this way is not 
rigorous and is still unresolved (Peacock & Heavens 1990). However the slight difference 
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in the results by adopting the different filters does not cause any essential change in the 
main conclusions drawn in the later sections. 

3.2 The effects of cooling and UV background 

The mass functions with the threshold value for v given in eq.(7) apply only to objects 
which can collapse and virialize before a given redshift. In order to form stars, however, 
the hot interstellar gas with virial temperature, T ~ 10 5 ~ Q K, must cool via radiative 
processes. Furthermore, in the presence of UV background, the gas must be dense enough 
to form a neutral core. 

These two requirements on the virialized objects can be expressed in terms of the 
mass-dependent thresholds, v cool (M), below which the objects cannot cool before a given 
z in the absence of UV background, and v uv (M), below which they cannot form self- 
shielded cores against UV photons. We define v cool (M) for the given total mass M, such 
that the cooling time given in eq.(l) (without the effect of photoionization) satisfies the 
following condition, 

tcool £ -|-[(1 + Z)- 3 / 2 - (1 + Z C )-" 2 1 (12) 

where z c is the redshift of collapse. At the high mass end for M ^ 1O 1O M0, line cooling 
by metals produced in galaxies formed earlier is important. To include the effect of metal 
cooling, we have adopted the simplification, as in Peacock and Heavens (1990), that for a 
cooling function dominated by recombination below T ~ 10 8 K, the collapsed redshift can 
be scaled as 

(1 + z c ) = (1 + z)(l + MjM cool fl\ (13) 

where 

M cool ~ 1.68 x 10 13 / m M @ . (14) 

The parameter f m denotes the metal content: f m — 1 for solar abundances and f m ~ 0.03 
for no metals. We have assumed f m — 1 for metals and 10% for the gas fraction. M coo i 
then is the typical mass scale below which cooling is important. Substituting the above 
expression for (1 + z c ) in eq.(7), the threshold for the density fluctuation in eq.(7) will have 
an extra factor for <5 C , i.e. 6 C — > 6 C (1 + M / M coo {) 2 ^ , indicating the increase of the threshold 
due to the cooling condition. At the low mass end for M ^ 1O 1O M0, since the effect of 
metalicity is trivial, we have determined v cool (M) based upon the zero-metal cooling rate. 
As we will show below, the condition to form a shielded core represented as dash-dotted 
curves in Fig.l is more stringent than the cooling condition. Our use of two different 
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cooling functions in two different mass ranges does not affect the final result for the mass 
function of galaxies. 

The effect of UV background on the threshold, v uv (M) was obtained in §2. Fig. 2 
shows that the shielded core mass (including both dark halos and baryons) is almost equal 
to the total mass for v(M) > v uv (M), and is zero for lower values of v(M). 

The new threshold value for density fluctuations which can self-shield and cool, is then 
given as 

v th {M) = max{v, v coo \ v uv ). (15) 

Fig. 4 shows the threshold values for collapse and cooling of objects by the present epoch 
z — 0, u uv for self-shielding in the UV flux J_2i = 0.1, and the adopted v t h given in 
eq.(15), as a function of core mass Mi. As discussed in §2, u uv is calculated on the 
assumption that the baryonic gas with the fraction, F — 0.1, in mass, collapses by a 
factor F -1 / 3 within a uniformly distributed dark halo. At both high-mass Mi £ 1O 12 M 
and low-mass Mi ;$ 1O 7 M0 ends, v t h is very large because of the cooling constraint. It 
is also evident that the condition to form a shielded core against UV photons essentially 
determines the v t h in the mass range 10 8 ^ Mi ^ 1O 12 M0, i.e. from dwarf to bright 
galaxies. This is clearly seen in Fig.4b where the dependence of v t h on the UV flux J_2i 
is shown at a fixed mass scale Mi — 10 11 Mq. It shows that v t h increases with J_2i- 

Having obtained v t h(M) for several constraints, the number density of objects as a 
function of core mass Mi can be calculated based upon eq.(9) and eq.(ll) for the general- 
ized PS and peaks formalisms, respectively. 

3.3 Mass function at z—0 

In Fig. 5, we show the mass functions of baryon masses Mi at z — under the above 
conditions, based upon the PS and peaks formalisms. The UV flux, J_2i, is assumed 
to be 0.1. The mass functions for collapsed objects or peaks have the power-law tail 
dn/dMi oc M^ 2 on small mass scales. The exponential cutoffs at high-mass ends exist 
simply because the objects with 5 < S c for all larger radii cannot collapse by the present 
time. The mass function in the peaks formalism has a slightly steeper slope at low- 
mass ends than that in the PS formalism because of counting the small-scale peaks which 
are actually within larger systems. The amplitude of mass function at low-mass ends is 
obviously related to the efficiency of merging. 

The mass function with cooling condition denotes the objects which have cooled by the 
present time. In both formalisms, the numbers of high mass objects are drastically cut off 
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due to ineffective cooling at these mass scales. In the PS formalism, the number density at 
intermediate mass Mb ~ 1O 12 M is increased, because the ineffective cooling of more mas- 
sive objects inhibits the merging of objects at this mass scale and their survival increases 
their number density (see Peacock & Heavens 1990). For low-mass scales Mb ?S 1O 7 M , 
where the virial temperature is below T ~ 10 4 K and cooling is negligible, there appears 
a cutoff in both formalisms. The ineffective cooling in these mass scales leads to the high 
threshold value v (Fig.4a) for density fluctuation, resulting in the suppression of low-mass 
galaxy formation. According to J. Peacock (private communication), this decrease of the 
probability in the PS formalism that low-mass objects can cool is compensated by the high 
probability for larger objects, and this is the origin of the bump just above this cutoff of 
the mass function — as also happens just below the high-mass cutoff. Note that even if the 
present recipe fails to derive mass function at low-mass ends, such very low-mass objects 
are observationally unimportant due to the selection effects (Yoshii 1993). 

Under the cooling condition, there is a slight dip in the PS mass function at Mb ~ 
10 s ' 5 Mq. This is artificially produced in our method, because our v t h given in eq.(15) and 
Fig. 4a has a discontinuous change from the cooling-condition dominant (low-mass) scales 
to the collapse-condition dominant scales (at Mb ~ IO^Mq). Probably, nature would 
avoid such a discontinuity in e.g. due to the effects of non-uniform and non-spherical 
density distribution of gas, and smooth out the mass function. However such a feature in 
the mass function does not affect the main conclusions below. 

Imposing the more stringent condition of forming self-shielded cores in the presence 
of UV background leads to further decrease in the mass function. Table 1 and 2 give the 
number density as a function of the UV flux for a few mass scales. It is evident that at 
the high-mass end, the number of objects which can form a core is smaller than that of 
objects which can cool by the present time. This is due to the efficient photoionization of 
the low-density gas in the massive objects (see Fig. 4). For masses M b £ 1O 13 M , however, 
the cooling condition becomes more important, giving rise to a sharp cut off. At low-mass 
end, the peaks formalism predicts a decrease in the number density, whereas it does not 
change much in the PS formalism. As a result, the slope of the mass function in PS 
formalism becomes steep due to the UV background, while it is slightly steepened in the 
peaks formalism. An artificial irregularity in the PS mass function, as happened in the 
case without UV, also appears at Mb ~ 1O 8 M0 in Fig. 5a, where our v t h given in eq.(15) 
and Fig. 4a has a discontinuous change from the cooling-condition dominant (low-mass) 
scales to the shielding-condition dominant scales. 
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To estimate quantitatively the change of slope, we assume a power-law form, 
dn/dlnM b oc for the mass function at low-mass ends 1O 8 M ^ Mb ^, 1O 1O M , and 

fit the derived mass function, a is a constant index denoting the slope of mass function at 
low-mass ends. In table 3, we show the values of a as a function of J_2i- ot is apparently 
a decreasing function with J_2i- It is interesting to note that the resultant shape and 
slope of tions in the UV background are similar for both formalisms, albeit with higher 
amplitude in PS-formalism. 

3.4 Evolution of the mass function 

In order to calculate the evolution of mass function in the presence of the UV back- 
ground, we need to express the possible evolution of the UV flux, J_2i , as a function of z. 
As we mentioned in §2, several observations and calculations suggest that J_2i is almost 
constant before z ~ 2, and rapidly declines at lower redshifts. Since the knowledge of J_2i 
at intermediate redshifts is still uncertain, we have adopted a simple functional form for 
J_2i(z): 

T f ^ \ J-2i(z = 2) = const. for z > 2 

J - 2l(z) -\j_ 21 (z = 2)(^) a forz<2. (16) 

We adopt 3 characteristic cases for the evolution of J_2i , defined by ( J_2i (z — 2), J-21 (z — 
0)) = (1.0,0.1), (3.0,0.1), and (3.0,0.05). The parameters a's in eq.(16) are then equal to 
2.10, 3.10, and 3.73 respectively. 

In Figs. 6 and 7, we show the redshift dependence of mass functions in the range 
< z < 3 for PS and peaks formalisms, respectively. The collapse factor of baryonic 
gas is again i* 1-1 / 3 for the cases with UV background. The major difference in the mass 
functions between the two formalisms is that when there are no constraints due to cooling 
and UV background, the amplitude of mass function in PS formalism is a decreasing 
function with time (Fig. 6a), whereas the number of collapsed peaks is always increasing 
(Fig. 7a). This difference is due to the explicit consideration of merging in the former case. 
The characteristic mass M*, (where the power-law mass function turns to an exponential 
fall-off at high-mass end), increases in both formalisms, since larger systems collapse at 
later times. 

However, the effect of photoionization tends to smear out such a difference (Figs. 6b, 
6c, 6d, 7b, 7c and 7d). The evolution of mass function is greatly influenced by UV, 
especially in the mass range 10 8 < M b < 1O 1O M : the amplitude of mass function at a 
given mass scale tends to increase with time after z — 2, even in PS formalism. The 
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slope is slightly steeper at higher z, because the lower- density gas in more massive objects 
is more vulnerable to photoionization. Before z — 2, the shape of mass function does 
not change with time owing to the constancy of the UV flux at these epochs. The rate of 
increase in the amplitude is dictated by the rate of decrease in UV, i.e. the evolution of the 
mass function in time is greater for larger J-21 at z — 2 (for the same value J-21 (z — 0)). 

The effects of UV upon the mass function obtained here do not depend upon an im- 
posed bias parameter b. If we decrease the value of b from 2 to 1, the threshold value of 
density fluctuation, v t h, is decreased as shown in Fig. 4(b), and the amplitude of mass func- 
tion is accordingly increased. However, the steepened slope and the decreased amplitude 
of the mass function due to UV background are unchanged by this parameter. 

4. DISCUSSION 

4.1 Galaxy luminosity function 

To compare our results for the mass function with observations, one needs to derive the 
luminosity function from the mass function. The mass-to-luminosity ratio M/L, however, 
is generally not a constant in mass and time and detailed modeling of star formation and 
gas dynamics such as gas loss from galaxies are required to derive the actual luminosity 
function (e.g. Lacey & Silk 1991; White & Frenk 1991). Such a derivation of the luminosity 
function is therefore very model dependent. 

The possible model-independent effects of the UV background upon the luminosity 
function are: (i) the number density of galaxies decreases due to UV, if one adopts the same 
star-formation law for both models with and without UV; (ii) the slope of the function 
would change at a mass scale, say M s , where the UV-dominant regime for Mf, < M s turns 
to the cooling-dominant one for Mf, £ M s (M s ~ 1O 12 5 M0 in Fig. 5). The mass-dependent 
M I L would not change this property, because it is unlikely that this ratio varies by an 
order of magnitude in the narrow range of mass scales near M s , as is the case for the mass 
function (note the log-scale ordinate in Fig. 5). We will further discuss the effect (i) later 
in conjunction with the galaxy number counts. 

4.2 A biasing mechanism 

The necessity of some physical biasing mechanism in galaxy formation has been argued 
by various authors in trying to reconcile the CDM model with large-scale galaxy distribu- 
tion and velocity distribution. Rees (1984) had argued that photoionization should make 
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it 'difficult' for galaxies, specially those with v ^ 1, to collapse. We have been able to 
quantify the 'difficulty' in terms of the threshold amplitude of fluctuations, v t h(M). As 
Fig. 4b shows, the threshold amplitude is indeed increasing with J-21, which results in the 
suppression of the number density of baryonic cores (Table 1). We therefore have a natural 
biasing of the galaxies relative to the dark halos. For instance, according to Fig. 4b at the 
fixed mass scale Mb — 10 11 Mq, the value of Vth{z — 0) for J-21 — with a bias parameter 
b — 2 is comparable to that for J-21 — 0.05 without artificial bias, b — 1. Furthermore 
as is shown in Fig. 4a, the threshold amplitude tends to increase more at higher masses 
due to the UV background, simply because the low density gas at high mass scale is easily 
photoionized. 

Therefore, the presence of the UV flux in the background provides a biasing mecha- 
nism for the formation of galaxies: the effect is more at brighter galaxies, enhancing the 
clustering. This results in the preferential formation of bright galaxies in clusters and 
superclusters. It is also possible that the spatial distribution of quasars, which are the 
main contributors to the UV flux, introduces the modulations in the galaxy distribution 
(Efstathiou 1992). Bahcall and Chokshi (1991) suggested that quasars are preferentially 
located in small groups of galaxies (JV = 10 ~ 30), while according to Mo and Fang (1993), 
quasars trace the same large-scale structure of galaxies and cluster of galaxies. If quasars 
are located in high-density environment of galaxies, the UV flux there would be stronger 
than that in the field. This will thus result in steeper slope of the mass function in the 
richer environments due to stronger UV flux (see Table 3). The observed correlation be- 
tween the slope of luminosity function and the richness of group (e.g. Ferguson 1992 for a 
review) may be explained by such an inhomogeneity of the UV flux. 

It is true that our calculations for uniform spheres may not be realistic in the sense 
that one expects the density profile in the halos to be more like isothermal. The collapsing 
factor would be F~ l in that case, instead of _F -1 / 3 that we used. It is however not clear, 
without a detail computation that includes the density profile, how the density profile 
affects the results. 

4.3 Galaxy number counts 

The presence of the UV background has further implications in the interpretation of 
galaxy number counts. The number count at an apparent magnitude rri\ at z (where the 
subscript A denotes the wavelength dependency), with an absolute magnitude M\, is given 
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by the product of the luminosity function (f)(M\, z) and the volume element dV/dz of the 
comoving space: 

N(m x ,z) = ^(M x ,z)^, (17) 

where lo is the angular area in units of steradians over which galaxies are counted. Inte- 
grating over z, we obtain the differential number count of galaxies N(m\) as a function of 
m\. N(m\) is strongly dependent on the evolution of galaxies through M\(m\) relation 
and the evolving </>, and on the cosmological parameters through volume effects (e.g Yoshii 
& Takahara 1988; Rocca-Volmerange & Guiderdoni 1990). In order to explain the large 
number of blue faint galaxies revealed by the observation (Tyson 1988), several attempts 
have been made, which include: a merging-driven evolution with Q — 1 based upon the 
Schechter-form luminosity function (Rocca-Volmerange & Guiderdoni 1990); the effect of a 
large cosmological volume element (low Q and non-zero cosmological constant A: Fukugita 
et al. 1990); the bursts of star formation at a recent epoch (Lacey & Silk 1991; Babul & 
Rees 1992). The interpretation of galaxy number counts is thus highly dependent upon 
luminosity and number evolutions as well as cosmological parameters. 

The mass function based upon the hierarchical clustering of CDM predicts the steep 
mass spectrum dn/dMf, oc M^ 2 on small mass scales (Table 3; Lacey et aZ.1993), in contrast 
with the observed slopes in the field, a ~ —0.95 (Loveday et aZ.1992) and in the cluster of 
galaxies, a ~ —1.3 (Sandage et a/.1985). This suggests the excess of low-mass galaxies to 
explain the faint galaxy count, but contradicts with the present galaxy luminosity function, 
even considering supernova- driven mass-loss processes (Lacey et a/.1993). 

In the presence of the UV background, the present work suggests that the number 
density of objects which can collapse and cool by a given epoch is smaller than that 
without UV due to the photoionization effect. Although the accurate incorporation of 
star formation and gas ejection via supernova-driven galactic wind, which make the slope 
of luminosity function shallower, is beyond the scope of this paper, the effect of the UV 
background may decrease the amplitude of luminosity function by reducing the number 
of the pre-galactic objects. Furthermore, because of the decreasing UV flux with time for 
^ z ^, 2, the time-evolution of the mass function is not constant nor decreasing with time 
as is usually assumed in the models of number counts (Yoshii & Takahara 1988; Rocca- 
Volmerange & Guiderdoni 1990). The effect of photoionization is thus suppressing the 
number of galaxies, preferentially at earlier epochs (Efstathiou 1992). 
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Although we have investigated only the case of Q — 1 universe, these effects of the 
UV background on the mass function are essentially the same for other cosmological pa- 
rameters. It was proposed, using the Schechter luminosity function dn/dL oc L~ l on low 
luminosity scales L, that the strong number evolution of galaxies via merging, i.e. larger 
number of building blocks (premergers) at higher redshifts, can save the Q — 1 universe 
(Rocca-Volmerange & Guiderdoni 1990). However, our result of increasing mass function 
amplitude with time by the UV background implies a larger discrepancy between observed 
and theoretically predicted number counts: the number of premergers is smaller at earlier 
epochs. This opposite sense of time-evolution in mass function may not support the as- 
sumption of the strong decrease in the number of galaxies adopted in the number counts 
model. Note here that the presence of a bump in the mass function at low-mass scale 
Mb ~ 10 7 Mq in the PS formalism is not important in this discussion because such low- 
mass galaxies at high z would not be detected due to the selection effects (Yoshii 1993). 
Therefore, it is likely that the high density Q — 1 universe does not reconcile with the 
large number of faint counts. 

We note from fig. 5 that the decrease in the number density of the galaxies is small 
for reasonable values of J_2i- This means that invoking the decrease in J_2i after redshift 
of z ~ 2 to explain the emergence of a number of small mass galaxies, as was done by 
Babul and Rees (1992) to explain the galaxy counts at faint magntitudes, can be fraught 
with difficulties. Our results do not support the idea that only the decrease in the UV 
background radiation after z ~ 2 could trigger an epoch of starbursts at z ~ 1 in dwarf 
galaxies, as is required in the mechanism of Babul and Rees (1992). 

5. CONCLUSIONS 

We have investigated the effects of the UV background radiation at redshifts z — ~ 3 
on the mass function of galaxies and its evolution in CDM model. Imposing the condition 
that, in order to cool for star formation, objects should be able to self-shield against 
photoionization, we have determined the threshold amplitude of density fluctuations. We 
have used the Press- Schechter and peaks formalisms to calculate the abundances of objects 
and evolution of the mass function in time. 

We found that the number density of objects is reduced in the presence of the UV 
photons. For example, at z — 0, and for J_2i(z = 0) = 0.1, the number density of 
M}, ~ 10 11 Mq galaxies is less by a factor of ~ 2 compared to the case of no UV radiation. 
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Our results, however, do not support the claim that an epoch of starbursts at z ~ 1 is 
triggered by the decrease in the UV background radiation after z ~ 2, as has been argued 
to explain faint blue galaxies. We also found that the number density of low-mass objects 
increases in time, in both Press-Schechter and peaks formalism, due to the evolution of 
the UV flux. This means that the difference between the cases with and without the UV 
background is larger at higher redshifts (0 < z < 2). This result is contrary to the usual 
view voiced in the literature, in the context of merging driven models of galaxy evolution, 
that number density of small mass objects decreases in time. 

We have discussed the implications of our result with respect to galaxy number counts 
and the luminosity function of the galaxies. We found the effect of the UV background 
radiation to be large, and it seems important to include its effect in any realistic model of 
galaxy evolution. 



We are indebted to Drs. Yuzuru Yoshii, John Peacock, Alan Heavens, and Peter 
Biermann for their useful comments and discussions. We are also grateful to the referee 
for his excellent and constructive suggestions for the improvement of this paper. MC 
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TABLE 1 

Changes in the number density (Mpc -3 ) (PS formalism) 



M b (M Q ) J_ 21 = 0.0 OT 05 1.0 

10 9 2.03-10- 1 2.26-10- 1 2.29-10- 1 2.17-10- 1 

10 10 3.11-10- 2 2.81-10- 2 2.58-10- 2 2.26-10- 2 

10 11 4.56-10- 3 2.99-10- 3 2.22-10- 3 1.80-10- 3 

10 12 8.29-10 -4 2.87-10 -4 1.35-10 -4 7.71-10 -5 



TABLE 2 

Changes in the number density (Mpc -3 ) (Peaks formalism) 

M b {M Q ) J_ 21 = 0.0 0.1 0.5 1.0 

10 9 2.94-10 -1 2.00-10 -1 1.47-10 -1 1.21-10 -1 

10 10 2.71-10- 2 1.68-10- 2 1.09-10- 2 8.84-10- 3 

10 11 2.14-10- 3 1.24-10- 3 7.12-10- 4 5.50-10- 4 

10 12 l.OMO- 4 8.12-10- 5 3.60-10- 5 2.24-1Q- 5 



TABLE 3 

Changes in the slope of mass function at low-mass scales 10 8 < M b < 1O 1O M0 

Formalism J_ 2i = 0.0 OT 05 O 

PS -1.89 -1.93 -1.96 -1.99 

peaks -2.02 -2.07 -2.10 -2.12 
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Figure Captions 

Figure 1. Gas number density versus virial temperature or virial velocity. Thick solid 
curves: the cooling curves where tff — t coo i, for J-21 — 0, 0.1, 1. Dotted curves: the equi- 
librium temperature T eq where cooling rate is equal to heating rate due to photoionization, 
for J-21 — 0.1 and 1. Dash-dotted curves: the self-shielding curves against UV photons 
for J-21 — 0.1 and 1. Solid curves: the loci of density perturbation for la in the virial 
equilibrium normalized to unity at a comoving radius lS/i^Mpc, with a bias parameter 
b — 1. The corresponding parallel solid curve (separated by an arrow) denotes the gas 
cloud after the density increase by a factor 10 due to a dissipative contraction. Dashed 
lines: the total-mass constant lines, assuming a baryon fraction of 0.1. 

Figure 2. Mass of a self-shielded core against UV photons, M c , as a function of total mass 
Mt (including dark matter for a baryon fraction of 0.1) for several values of the threshold 
is. Solid lines: for J-21 — 0.1. Dotted lines: for J-21 — 1.0. Attached numbers denote the 
value of v. 

Figure 3. Rms density fluctuation A as a function of the top-hat scale R for two filters, 
top-hat filter (solid line) and Gaussian filter (dotted line). Both curves are normalized so 
that A is unity at the top-hat lengthscale lGh^Mpc with h$Q — 1. 

Figure 4. Threshold of density fluctuation, is, affected by the UV background, (a) The 
mass-dependency of is at z — with b — 2, for several constraints: dash-dotted line: is for 
collapse before z — 0; dotted line: is for cooling before z — (without UV); dashed line: is 
for forming a self-shielded core against UV photons represented by J-21 — 0.1; solid line: 
adopted is to estimate the mass function, (b) The dependence of is on the UV flux J-21 
at the mass scale Mj, = 10 11 Mq, with bias parameters b — 1 and 2. Solid lines: for z — 0. 
Dotted lines: for z — 2. 

Figure 5. Number density of objects as a function of baryonic mass Mf, at z — 0, with a 
bias b — 2, in the Press-Schechter formalism (a) and the peaks formalism (b). Solid lines: 
mass function for collapsed objects at the present epoch. Dotted lines: mass function 
for cooled objects at the present epoch. Dashed lines: mass function for cooled and self- 
shielded cores against UV with J-21 — 0.1. 

Figure 6. Time-evolution of the mass function predicted in the Press-Schechter formalism 
with b — 2, at z — 3 (dash-dotted lines), z — 2 (dashed lines), z — 1 (dotted lines), and z — 
(solid lines), (a) mass function for all collapsed objects, (b) the case for J '-2i(z — 2) — 1.0 
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and J_ 2 i(z = 0) = 0.1. (c) the case for J- 2 i(z = 2) = 3.0 and J- 2 i(z = 0) = 0.1. (d) 
the case for J_2i(z = 2) = 3.0 and J_2i(z = 0) = 0.05. The adopted evolution of J-21 is 
given in eq.(15). 

Figure 7. The same as Fig. 6 but for the peaks formalism. 
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